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Context. The dipole modes are the most common geometry of oscillations in roAp stars inferred from photometric 
measurements and are therefore of special interest for asteroseismic purposes. 

Aims. We present a theoretical and analytical study of the light curves associated with dipole {£ — 1) pulsations of 
roAp stars in the framework of the revisited oblique pulsator model. 

Methods. We describe the light curves in terms of the inclination and polarization of the elliptical displacement vector 
of the dipole modes. We study the influence of the magnetic field and rotation on the shape of these light curves for 
both amplitudes and phases. 

Results. Despite the inclination of dipole mode with respect to the magnetic axis, we find that the dipole mode can have 
maxima that are in phase with the magnetic maxima. We apply our formalism to the well-known roAp star HR3831 
(HD 83368) to derive its mode properties. Our results are similar to those obtained by time-series spectroscopy. We also 
consider the cases of three other roAp stars, HD6532, HD 99563, and HD 128898 (aCir). 

Conclusions. We demonstrate that the formalism of the revisited oblique pulsator model is adequate to explain the 
properties of the photometric light curves associated with dipole modes in roAp stars. In addition, we show that the 
coincidence of pulsation and magnetic extrema can also occur for inclined modes with respect to the magnetic axis. 
With the stars considered in this paper, we conclude that the polarization of the modes present in roAp stars are quasi 
linearly polarized. 
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1. Introduction 

It is well-known that the rapid light variations in Ap stars 
are caused by trapped magneto-acoustic modes. The first 
detection of these oscillations was made by Kurtz (1978, 
1982), who clearly showed the presence of high radial order 
modes (n sa 20 — 30) with periods comparable to those of 
the solar-like stars; these periods are now known to range 
from 6 to 21 min. The photometric amplitudes of the roAp 
stars, with values of about 0.1 — Bmmag, are vastly higher 
than those of the solar-like oscillators. It is now common to 
refer to these objects as roAp stars, for rapidly oscillating 
Ap stars. To date, over 40 roAp have been detected, Kurtz 
et al. (2006) hsting 40 of them. 

We can distinguish three fundamental properties of the 
photometric variability of these stars: 

1. The light curves of many roAp stars show a double 
modulation. The rapid one, with a period of a few 
minutes, arises from the pulsation modes. The long 
one, with a period between days and decades, is the 
rotation period of the star. To explain the presence of 
the long-term variation, Kurtz (1982) proposed that 
the pulsation mode axis is aligned with the magnetic 
axis, which is itself inclined to the rotation axis so that 
the observer sees the pulsation modes from an aspect 
that varies with rotation. This is the oblique pulsator 
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model. 

2. The second main property of the roAp stars is that 
the maxima of the rotational amplitude modulation 
of the pulsations occur at the same rotation phases 
as that of the magnetic maxima. This observational 
finding led Kurtz (1982) to justify the inclination of the 
modes with the rotation axis, simply by considering 
a pulsation axis aligned with the magnetic one. In 
the magnetic Ap stars, it is indeed well-known that 
the magnetic and rotation axes are inclined with 
respect to each other (e.g. Landstreet & Mathys 
2000); this is the oblique rotator model (Stibbs 1950). 
Kurtz's explanation looked reasonable given the huge, 
kG-strength magnetic fields observed at the surface 
of these stars, at least two or three orders of mag- 
nitudes higher than the global magnetic field of the Sun. 

3. Finally, the last main property derived from photom- 
etry is that for some of the roAp stars - those with 
geometries such that both magnetic poles come into 
view over the rotation period - the light curves show 
a phase jump of ±7r radians as the pulsation amplitude 
passes through a minimum. This property is most easily 
explained when the geometrical structure of the mode 
consists of two hemispheres with opposite phases, say 
a dipole {£ = 1) mode, although other, higher spheri- 
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cal degree modes may show the same effect for certain 
viewing geometries. 

Over the past two decades, many theoretical investiga- 
tions have studied the properties of the roAp stars. The 
purely geometrical model of Kurtz was improved by con- 
sidering the effects of both the magnetic field and rotation 
(Dziembowski & Goode 1985; Kurtz & Shibahashi 1986; 
Shibahashi & Takata 1993; Takata & Shibahashi 1995). 
However, despite these improvements, the magnetic field 
treatment in all these studies was oversimplified. They in- 
deed considered a perturbative treatment to solve the mag- 
netic effects in the problem, which is unrealistic because in 
the outer layers of the star the magnetic pressure is higher 
than the gas pressure. 

In these models, the modes were found to be aligned 
(or nearly aligned) with the magnetic field, in agreement 
with Kurtz's original proposition of aligned pulsation and 
magnetic axes. The reason behind this alignment is that 
for high overtone pulsations, the magnetic perturbation is 
much larger than the Coriolis one. 

On the basis of both photometric observations and the- 
oretical properties, the picture of the pulsations in roAp 
stars over the past two decades has been that of dipole os- 
cillations with the pulsation axis aligned with the magnetic 
axis. The inclination of these with respect to the rotation 
axis thus explains the observed long-term modulation of 
the pulsation amplitude. 

Bigot & Dziembowski (2002) proposed an alternative 
description of pulsations in roAp stars. They found that 
the axis of pulsation can be inclined by a large angle to the 
magnetic axis. This result comes from a degenerate per- 
turbative treatment of the rotation effects for which they 
added the quadratic effect of the centrifugal force in the 
problem. This second order effect of rotation was consid- 
ered negligible in previous studies. For high radial order 
modes, such as those in roAp stars, the centrifugal effects 
can be comparable to the effect of the magnetic field. These 
improvements are crucial for quantitative studies and can 
change the inclination of the mode. 

For relatively moderate fields, unlike in the previ- 
ous studies, there is no dominant perturbation acting 
on the mode, but rather a balance between magnetic 
and rotational perturbations. In their formalism. Bigot 
& Dziembowski (2002) showed that in the case of dipole 
modes the pulsation axis lies in a plane during the pulsation 
cycle, with a displacement vector moving with an elliptical 
motion. These dipole modes are fully characterized by three 
quantities: the inclination of the mode plane, the elliptic- 
ity of the displacement vector, and the frequency. Bigot and 
Dziembowski applied their model to the roAp star IIR3831 
and based on the magnetic measurements (Bagnulo et al. 
1999), they showed that, for this star, the mode axis is well 
inclined to both the magnetic and rotation axes. 

Considerable progress has since been made in the un- 
derstanding of pulsation in roAp stars by analysing spectro- 
scopic time-series observations. In particular, those lead to 
the possibility of probing the vertical structure of pulsation 
through the atmospheric layers by analysing the pulsation 
through different lines that form at different atmospheric 
depths. The multiplets observed in spectroscopy do not look 
like those derived by photometry, but have instead stronger 
non-dipolar components. Kochukhov (2006) derived a new 
geometrical aspect of the mode of IIR3831 using radial ve- 



locity measurements and found no significant inclination 
with respect to the magnetic axis. Therefore, he returned 
to the original proposition of Kurtz (1982) and thus ques- 
tioned the results of Bigot & Dziembowski (2002). 

The main objective of this paper is to provide an ob- 
servational counterpart to the model developed in Bigot & 
Dziembowski (2002). We derive a formalism to calculate 
analytical light curves associated with an arbitrary combi- 
nation of dipole modes (m = —1, 0, +1). The relatively sim- 
ple geometrical model that we propose is particularly use- 
ful for analysing photometric pulsation data. We see that 
it closely agrees with both photometric and radial velocity 
observations. 

In section 2, we give an analytical expression for the 
light curves of dipole modes in terms of the formalism pro- 
posed by Bigot & Dziembowski (2002), i.e. the mode incli- 
nation and polarization. In section 3, we discuss the cor- 
relation between the dipole mode inclination and coinci- 
dence of pulsation and magnetic extrema. In section 4, we 
consider the effects of a changing magnetic field on these 
light curves and how it affects their amplitudes and phases. 
Finally, in section 5 we apply these new considerations to 
the cases of four roAp stars IIR3831, HD6532, a Cir, and 
HD 99563. We see that our formalism agrees well with ob- 
servations. For the particular case of HR3831, in the light 
of a new measurement of the inclination angle of the ob- 
server, we show that our formalism leads to results that 
are consistent with those of Kochukhov (2006) within the 
formalism presented in the oblique pulsator model of Bigot 
& Dziembowski (2002). 

2. The light curve variations of dipole modes 

2.1. A general expression for non-radial oscillations 

The formalism for light variations in the case of non-radial 
oscillations is given in Dziembowski (1977) (see also Buta 
& Smith 1979). We follow the formulation of Dziembowski, 
but with an additional approximation specific to roAp 
stars. In the case of high radial orders, the main contri- 
bution of the fluctuations to the luminosity SL comes from 
the variations in the radiative flux 6J-. We can reasonably 
ignore the effects of any changes in the surface element, in 
terms of both area and orientation. Thus, the light varia- 
tions expressed at the stellar surface may be written as 

= ^HlJ')lJ'd^^d(j), (1) 

where ^ = cos and h{fi) is the so-called normalized limb 
darkening function (e.g. Mihalas 1978). 

The fluctuation SJ- is the unknown quantity that re- 
mains to be determined. In roAp stars, the modes are 
strongly influenced by both magnetic and inertial forces. 
The joint effect produced by these two perturbations is non- 
spherically symmetric, leading to a distortion of the sur- 
face amplitude of the modes. Mathematically, this means 
that (5J-" is written as a linear combination involving spher- 
ical harmonics of different degrees i and corresponding az- 
imuthal orders m = [—£, £]. The eigenmodes are solutions of 
an eigenvalue system. The eigenvalues are the frequencies, 
and the eigenvectors are the coefficients of the spherical 
harmonic expansion. 
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In this context, Kurtz (1992) proposed that the light 
variations of roAp stars could be interpreted in terms of ax- 
isymmetric modes (m = 0) about the magnetic field, whose 
angular dependence is described by a linear combination of 
i (= 0, 1,2,3). This combination of spherical harmonics was 
used to take into account the distortion produced by the 
magnetic field. Dziembowski & Goode (1996) and Bigot et 
al. (2000) showed in the case of dipole magnetic fields that 
the distortion of the eigenvectors is quite weak as long as 
the field strength remains < 1 kG and can be regarded as 
a small correction to the problem. Cunha & Gough (2000), 
Saio & Gautschy (2004), Saio (2005), and Cunha (2006) 
explored the effects of stronger fields and showed that the 
surface amplitudes of the modes are strongly modified by 
kilogauss field strengths as found in roAp stars. However, 
photometric observations show dominant triplets, which 
mean that either the magnetic distortion is small or the 
geometric cancellation (integration over the stellar disk) of 
components £ > 1 leads the observer to see the distorted 
mode as a dipole. In the following, we assume a negligible 
magnetic distortion in order to develop an analytical for- 
malism, keeping in mind that it might not be applicable to 
all roAp stars. 

In this respect, the perturbation of the radiative fiux 
may be written as 



(2) 



where am are the coefficients of the coupling and lo the 
frequency. These quantities are, respectively, the eigen- 
vectors and eigenvalues of an homogenous system (e.g. 
Dziembowski & Goode 1985; Kurtz & Shibahashi 1986; 
Bigot & Dziembowski 2002). They depend on the magnetic 
and rotation strengths and orientations. The quantity fi is 
a function that characterizes the radial dependence of 6J- 
and is expressed here at the stellar surface. The coefficient 
e is unknown in the linear treatment of the oscillations, but 
is a small parameter so the perturbation formalism is valid. 
The Y™ are the usual spherical harmonics expressed in the 
rotation reference system, i.e. fiR = 1 corresponds to the 
rotation axis. 

We emphasize that the present approach differs from 
that of Kurtz (1992) since we consider that the distortion 
of the mode is described by a linear combination, not for £, 
but for the azimuthal order m. We stress that because of 
the inclination between the magnetic and rotational pertur- 
bations the coupling of modes for different azimuthal orders 
m can be strong and cannot be ignored in roAp stars. 

To proceed further and compare with observations, we 
have to write the flux variations described by Eq. ^ in the 
inertial reference system whose polar axis coincides with 
the line-of-sight of the observer. This operation is done via 
the well-known relation for spherical harmonics 



(3) 



Moreover, the observer being in an inertial frame, we have 
to transform the coordinate system as — >■ (pR — VLt, where 
51 is the rotation frequency of the star. Therefore, after 
substituting Eq.[3]into Eq.[51 and integrating over the stellar 
disk (only j = remains in Eq.[3]) the expression for the 
fluctuation in luminosity is given by 
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Since we work at a given value of £ and consider only 
relative fluctuations in the following sections, we are not 
interested in the value of the constant in front of the sum 
in Eq.m and just consider 
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(7) 



This expression clearly shows the double modulation ob- 
served in roAp stars, a rapid one of period 27r/w due to 
the oscillation and a long one of period 2^/0. due to the 
rotation of the star. In the frequency domain, it leads to a 
(2£+ 1) multiplet of amplitudes Am equally spaced by the 
frequency of rotation. 

2.2. Light curve variations of dipole modes 

The dipole modes have the most common geometry ob- 
served in photometry in roAp stars. They are therefore of 
special interest to asteroseismic studies. In the rest of this 
work, we will consider the fluctuations of luminosity in Eq. 
[7] for generalized dipole (£ = 1) modes, following the for- 
malism developed in Bigot & Dziembowski (2002). We will 
present an analytical formulation of the light curves for 
both amplitudes and phases. 



2.2.1. Properties of dipole modes 

As shown in Bigot & Dziembowski (2002), the dipole mode 
axis lies in a fixed plane during the pulsation cycle with an 
elliptical displacement vector. The orientation and shape 
of the ellipse are defined by its two axes: the A"-axis in the 
plane defined by the vectors (B, fi) and the perpendicular 
y— axis. The inclination of the X-axis of the ellipse with 
respect to the rotation axis is given by 



7 = arctan 



a+i — a_i 



(8) 



where {iiltiPl) are the coordinates in the new (observer's) 
reference system, i the angle between the rotation axis and 

the line-of-sight, and d^^, coefficients that are expressed 
in terms of the Jacobi polynomials (e.g. Edmonds 1960). ^ 



The elliptical displacement of the mode is characterized by 
the polarization angle 



a+i - a-i 



■ sin 7 



(9) 
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Dipole axis 




Fig. 1. The oblique pulsator model for dipole modes (Bigot & Dziembowski 2002). Dmlng the pulsation cycle, the axis 
of the dipole mode lies in a fixed plane that is inclined to both the rotation and magnetic axes. The shape of the 
displacement vector is elliptical. The two axes of the ellipse are the X— axis in the (B, Jl) plane and the perpendicular 
y— axis. The x-axis lies in the rotation equatorial plane. The line-of-sight of the observer makes an angle a(t) with the 
former, and a{t) with the latter. The magnetic axis is inclined by the angle X(t) to the line-of-sight. 



which measures the ratio of the two axes of the ellipse. 
The projected amplitudes of the displacement vector onto 
the X- and y-axes are then cos'0 and sin'0, respectively. 
We note that 1^ = corresponds to a mode that is linearly 
polarized in the plane (B,r2), whereas ip = ±7r/2 corre- 
sponds to a mode that is linearly polarized along the y-axis. 

The special cases are: 

— (I7I = 0,ip = 0) for a mode linearly polarized along 
the rotation axis, i.e. an uin = mode in the rotation 
system; 

— (I7I = 7r/2, ip = ±7r/4) for a mode circularly polarized in 
a plane orthogonal to the rotation axis, i.e. an niR = ±1 
mode in the rotation system; 

— (I7I — — 0) for a mode linearly polarized along the 
magnetic axis, i.e. an ms ~ mode in the magnetic 
system; 

— (I7I = 7r/2 — I3,ip = ±7r/4) for a mode circularly polar- 
ized in a plane orthogonal to the magnetic axis, i.e. an 

= ±1 mode in the magnetic system. 



2.2.2. Amplitudes and phases of the light curves 

The light curve associated with a dipole mode may be ob- 
tained from Eq.[7]with £ = 1. We write it in a form similar 
to that in Kurtz (1992) 



SL 
T 



but with different coefficients A{t) and B{t), 



A{t) = oqCOSz -|- "'^^ " ^ sinicosf^t, (11) 
v2 

and 

B(t) = ^±^-ti^sinisinm. (12) 

These expressions can be written in terms of 7 and ip as 
A{t) = cosip (cos 7 cos i + sin 7 sin i cos ilt) , (13) 
and 

B{t) = sin ip sin ism ^It, (14) 
where we have normalized each coefficient a„i (which are 



a? -I- ab- 



solutions for a homogeneous system) by 

From the observer's point of view, the relevant quanti- 
ties are the orientations of the X— and y— axes with respect 
to the line-of-sight, rather than the orientation with respect 
to the rotation axis (z-axis) . It is then more natural to con- 
sider the angles made by these two axes to the line-of-sight 
instead of the angles 7 and i. 

To define these angles, we consider the three directions 
of the problem that are defined by the three following unit 
vectors: 



cx A{t) cos ujt + B{t) sin cot, 



(10) 



1. the direction of the inertial observer 
Hobs = (sin i cos $(i), sini sin$(t), cos i), 
where $(t) ~ fit is the rotational phase; 



(15) 
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2. the direction of the X-axis 
nx = (sin7,0,cos7); 

3. the direction of the y— axis 
ny = (0,1,0). 



(16) 



(17) 



The angle a{t) between the X-axis of the ellipse and the 
line-of-sight is then defined by 

cos a{t) = nx-Hobs ~ cos i cos 7 + sin7sinicos(f>(t), (18) 

and the angle a{t) between the y-axis of the ellipse and the 
linc-of-sight is 



cos a{t) = ny.nobs = sin tsin$(i). 



(19) 



The expressions for the coefficients A{t) and B{t) reduce to 
simple forms 



A(<I>) = COS'(/)COSQ!($) 

i3($) = sin'(/)C0scr(<i>). 



(20) 
(21) 



These quantities then have a simple physical explanation: 
they correspond to the lengths of the X— and y— axes of 
the ellipse (cos'0 and simp, respectively, in the mode plane) 
projected onto the linc-of-sight. The amplitude and phase 
of the signal related to the dipole modes are finally only 
functions of three quantities: the eccentricity of the ellipse 
(ip) and its orientation with respect to the observer (a, a). 

Following Kurtz (1992), we write the light curve in 
terms of its amplitude R{t) and phase ^'(t) 



5L 

T 



(t) cx R{t) cos(a;i- ^'(t)). 



with 



and 



*(i) = arctan 



B{t) 



A{t) 



(23) 



(24) 



The behaviour of the amplitude R and phase for different 
mode inclinations and polarizations is represented in Fig. [51 

3. Extrema of the rotationally-modulated 
amplitudes of dipole modes and their 
coincidence with magnetic maxima 

The number and nature of the extrema of the rotationally- 
modulated amplitudes depend on the aspect of the dipole 
mode seen by the observer. In general, within one period of 
rotation, $ = [0, 27r], we count four extrema located at 



$ = [tt] and $ = $ J = ± 1 



1 



(25) 



where r = tan'!/'/sin7 and ^ = cot 7 cot i. The presence of 
four extrema of the envelope, despite the dipolar nature of 
the mode, is a direct consequence of the definition of -/?($) 
introduced in Ea.E^ 



In some special cases (polarization and/or inclination), 
the observer might see only two extrema. This may occur 
for example if |^| > — 1|. In that case, the observer sees 
extrema only at $ = 0[7r]. The inclination of the observer 
relatively to the mode axis is such that the same hemisphere 
is always seen as the star rotates. An illustration is given 
in the appendix in the special case of '0 = modes. 

Another case for the presence of only two extrema is 
found when the displacement vector of the mode is always 
perpendicular to the (B,J1) plane, i.e. when \tp\ = tt/2. In 
that case, there are no extrema at $ = [tt] , but only at 
$ = $0 = ±7r/2. 

In the following, we focus our discussion of the extrema 
located at $ = 0[7r]. The objective is to determine under 
which conditions the pulsations show maxima in phase with 
the magnetic maxima. 

3.1. The extrema of the magnetic field 

Owing to its inclination with respect to the rotation axis, 
the magnetic field is also modulated as the star rotates. The 
modulation depends on the orientation of its axis with re- 
spect to the line-of-sight of the observer. This describes the 
standard, well-known oblique rotator model (Stibbs 1950). 
Here we present this model in our notation. 

The direction of the magnetic field axis with the rotation 
axis is defined by the unit vector 



ne = (sin/3, 0,cos/3), 



(26) 



where (3 is the angle between the two axes. From the ob- 
server's line-of-sight, the magnetic axis is inclined by an 
angle A defined by 



(22) cos A($) = nB.Hobs = cos/3cosi -I- sin/3sinicos$(i). (27) 



Hence, the magnetic modulation due to the rotation of the 
star is 



(B) cx cosA($) 



(28) 



This magnetic variation is a sinusoid with two maxima at 
$ = and $ = TT. These maxima can be positive or negative 
depending on the value of A at $ = or $ = tt. 

The condition that has to be met for a positive maxi- 
mum at $ = and a negative minimum at $ = tt is 

|A(0)| < ^ < |A(7r)| or equivalently \l3-i\ < | < \l3+i\-{29) 

The first inequality ensures that there is a positive maxi- 
mum at $ = 0, i.e. |A(0)| < f , and the second one ensures 
that there is a second negative maximum, i.e. |A(7r)| > |. 



3.2. The extrema of pulsation at $ = [n] 

We now consider the nature of the pulsation extrema at the 
same rotation phases as those of the magnetic extrema. The 
first requirement to have magnetic and pulsation maxima 
in phase is obviously that the rotational envelope shows 
extrema at $ = [tt] . This means that the mode must have 
a non-zero amplitude in the magnetic plane (B,r2), which 
requires that |'0| 7^ tt/2 [tt]. If this condition is satisfied, 
the dipole mode always has two extrema in phase with the 
magnetic extrema at $ = 0[7r], since the X— axis of the 
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rototionol phase rototionol phase 



Fig. 2. Shapes of the light curves influenced by different inclinations and polarizations of the dipole mode. In both panels, 
the inclination of the observer is i = 70°. (Left panel) We consider a given polarization of the mode '0 = 5°, and compute 

curves for various inclinations of the mode plane: 7 = 5° ( — ), 25° (• • •), 45° ( ), 65° ( ), 85° (— • • • — ). It is clear 

that for all inclinations but one (7 = 5°), the phase jumps tt radians exactly at the amplitude minima since the modes 
show the two hemispheres during the rotation period (i.e., this verifies the conditions given in Eq. IA.2p . The position 
of the minima at are shifted by a change in the inclination of the mode, as one can check in Eq.[25] (Right panel) 
Shapes of the light curves influenced by different polarizations of the dipole mode. We consider a given inclination of the 

mode 7 = 45°, and compute curves for various polarizations ip of the dipole axis: -0 = 0° ( — ), 20° (• • •), 40° ( ), 

60° (— • — ), 80° (— ). In both panels, we have represented two periods of rotation. 



ellipse is within the magnetic plane. The question is then to 
know whether these extrema are maxima or minima. The 
answer depends on both the polarization ■0 of the dipole 
mode and its orientation toward the observer. 

At the rotation phase <& = [tt] , the angles between the 
X-axis of the ellipse and the line-of-sight are 



q;(0) = 7 — i and a{'K) = 7 + i. 



(30) 



The j/-axis of the ellipse is perpendicular to the line-of-sight 
((7 = ±7r/2) so its projected component becomes zero. The 
condition under which i?($) shows a maximum at $ = is 
easily derived by a simple analysis of Eq.[531 We find that 



cos ol[t:)) 



,,2 cosQ;(0)(cosa(0) 

tanV') < ^r-r-iTT^ • 

^ ' 2cos2cr(7r/2) 

A similar expression is obtained at the phase $ = tt, 

,2 cos a(7r)(cos a(7r) — cos a(0)) 

(tanw) < n — -, — TTT . 

^ ^' 2cos2a(7r/2) 



(31) 



(32) 



The conditions Egsl^ and are restrictive, in the sense 
that for a given position of the observer relative to the mode 
axis (i.e. X-axis at $ = [tt]), there is upper limit to the 
strength of the polarization allowed. In other words, the 
displacement vector has to stay close to the (B,ll) plane 
in order to show a maximum in phase with the magnetic 
one. 

It is important to note that these relations put 
some constraints on the polarization of the dipole mode. 
However, they allow modes that are (well) inclined to the 
magnetic field and have maxima in phase with the magnetic 
ones. Until now, it has been common to assume that the co- 
incidence between the pulsation and magnetic maxima im- 
plies that modes are linearly polarized along the magnetic 



axis, say = and 7 = /3. Even though this case does lead 
to such a coincidence, it is certainly not the only possibility. 
The relations in Eg s I5T] and [5^ indeed do not exclude modes 
with circular polarizations, ^ — ±7r/4, and even modes that 
are more elongated along the y-axis (7r/4 < |-0| < 7r/2). 

This result is not obvious: modes that have the largest 
amplitude of their displacement vector in a direction per- 
pendicular to the magnetic plane (B, $7), i.e. modes such as 
(I ■01 > 7r/4), can show some maxima in phase with the mag- 
netic maxima. This apparent contradiction can be easily 
explained because the observer does not see the real ampli- 
tude of the displacement vector but rather the projection 
of the X- and y-axes of the ellipse onto the line-of-sight. 
What really matters is the amplitude of the projection of 
the X-axis relative to the amplitude of the projection of 
the 2/-axis. The mode shows a maximum in phase with the 
magnetic maximum if the projection of its displacement in 
the plane (B, $7) is larger than the one perpendicular to 
this plane. 

To illustrate this, we consider an example: a mode with 
■0 = 55°, 7 = 50°, and i = 10°. This mode is more elongated 
along the y-axis than the X-axis, since sin'0 > cos-i/j. The 
inequalities in equations Eqs [31] and |32] are fulfilled, and 
the two extrema at $ = and $ = tt are maxima. The 
projection of the X-axis onto the line-of-sight is a($ = 0) = 
0.43, whereas the projection of the y-axis is fe(<I> = 7r/2) = 
0.14. Hence, the mode shows a maximum at $ = because 
the projection of the displacement vector component in the 
plane (B, $1) is larger than its component perpendicular to 
this plane (y— direction). 

The two extrema have unequal amplitudes because the 
two poles of the mode have different inclinations (q;(0) ^ 
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Fig. 3. Representation of the rotational variation of the light curves in amplitude and phase for three dipole modes on 
the six rows, and for three regimes of magnetic field in the three columns, /x = —0.2, —1.0 , and —5.0, where fi measures 
the ratio of the magnetic to centrifugal shifts of the pulsation frequencies. The amplitudes are normalized so that the 
maximum is equal to unity. In all cases, we choose i = 70° and /3 = 12°. The ratio of the Coriolis to centrifugal shifts of 
frequencies is x = 0.01. See the text for a description of the properties of the three dipole modes. The values (in degrees) 
of the inclination and polarization are given for each mode in the brackets (7, ip) above the panels. 



a(7r)) with respect to the line-of- sight. The ratio of their 
amplitudes is given by 



R{0) 



cos a{n) 



cosa(O) 



cos(i + 7) 



cos(i — 7) 



(33) 



This ratio is independent of i/j since at $ = [tt] the small 
axis of the ellipse is perpendicular to the line-of-sight (null 
projection). We note that for some special values of the 
inclination of the observer, i = 0[7r/4], this ratio equals 
one, regardless of the inclination of the mode. These cases 
are unfortunate for asteroseismology because they cannot 
provide information on the mode. 



4. Magnetic versus rotation effects on the light 
curves 

We now consider the influence of a varying magnetic field 
and rotation on the aspect of the light curves. At a given ro- 
tation frequency f2, a change in the magnetic field strength 



modifies both the inclination and polarization of the mode, 
which thereby changes its aspect seen by the observer. The 
results on the combined effects of a magnetic field and rota- 
tion on the dipole mode property are discussed in detail in 
Bigot & Dziembowski (2002). To relate the present work to 
these results, we keep the same values for the mode parame- 
ters and the same notation as Bigot & Dziembowski (2002). 
The magnetic field strength is parameterized by the value 
of fi, which measures the ratio of the magnetic to centrifu- 
gal shifts in frequencies. We emphasize that no assumption 
is made regarding the method used to solve the magnetic 
shift of frequencies that enters into the /i parameter. The 
values of > 1 correspond to a predominant magnetic 
field regime, and \fj,\ < 1 corresponds to a regime where 
the centrifugal force dominates. In Fig.|31 we use /z = —0.2, 
— 1.0 and —5.0. The ratio of Coriolis to centrifugal shifts of 
frequencies is kept to a constant value of 0.01, which is close 
to the real value corresponding to the roAp star HR3831. 

For dipole modes, the eigenvalue system leads to three 
orthogonal solutions that are represented on each row of 
Fig. 121 The aspect of these three modes is strongly modified 
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when the magnetic field strength increases, as one can see 
from left to right : 

1. In the low magnetic field regime {\fj,\ < 1, left column), 
each eigenmode is almost described by a single spheri- 
cal harmonic in the rotation system (mji — 0, ±1). The 
mode that is almost an mn = is not modulated be- 
cause its axis nearly coincides with the rotation axis. 
The phase of the modes that are almost = ±1 
modes show this characteristic straight line variation 
for the phase. 

2. The most complicated situation arises when magnetic 
and rotation effects are comparable (|/i| = 1, middle 
column), where the coupling between m-components of 
the mode is the strongest. In that case, it is impossible 
to identify a mode with a single value of m. 

3. In the case of strong magnetic field (|^| > 1, right col- 
umn), the modes are again almost described by a single 
value of m but now in the magnetic reference system. 
The one that is almost described by an niB = mode 
shows a constant phase except for phase reversals at 
amplitude minima. This jump corresponds to the mo- 
ment when the node of the oscillations passes through 
the plane (i, ft). The two other orthogonal modes are al- 
most circularly polarized (ip = ±7r/4) and tend to show 
characteristic monotonic linearly increasing or decreas- 
ing phases. 

For a given mode (a given row in Fig. [3]) , the light curve 
aspect changes drastically with the value of /z. This is be- 
cause the mode tends to be aligned with {uib = 0), or 
orthogonal (|tob| = 1) to, the magnetic axis as the field 
strength increases. This changes both the numbers and the 
amplitudes of the maxima of the rotational envelope, as well 
as the polarization. The mode that is almost an \mfi\ = 1 
mode in the rotation frame for fi < I becomes an = 
mode in the magnetic frame when /i > 1. Its phase curve 
changes from nearly straight lines to constant phase with 
jumps. 

We emphasize that there is a saturation effect once the 
magnetic field (equivalently /x) reaches a certain value 
4 — 5). Once the field is large enough, the three eigenmodes 
are either aligned or perpendicular to the magnetic axis 
and an increase in the field strength no longer has an effect 
on the aspect of the mode seen by the observer. Therefore, 
it becomes impossible to constrain the field strength above 
a certain limit by fitting the relative mode amplitudes. A 
tighter constraint is provided by the frequency spacing. 

5. Applications to four roAp stars 

5.1. Application to HR3831 

We apply these considerations to the roAp star IIR3831 
for which we have good time measurements for both mag- 
netic and pulsation data. This star pulsates in a single mode 
with a frequency of a;/27r = 1427/iHz (Kurtz ct al. 1997a). 
From photometric data, the mode appears essentially as a 
dipole since there is a clear phase reversal of n at the min- 
ima of the rotational envelope. Equivalently, this appears 
in the frequency domain as a dominant triplet with peaks 
that are exactly separated by the rotation frequency, and 
is clear evidence of a dipole mode pulsation. Extra £ ^ 1 
components are present in the spectrum, which indicates 



a distortion of the dipole by the magnetic field (Kurtz et 
al. 1997a), but they appear small in photometry, possibly 
because of either a slight influence of the magnetic field or 
a geometrical cancellation (Saio & Gautschy 2004 sec. 5.4) 
of the distorted mode. 

The radial velocity measurements of Kocliukhov (2006) 
indicate that the Nd ill lines have a strong £ = 3 component, 
much larger than in photometry, which is in strong disagree- 
ment with the mode geometry assumption of £ = 1 made in 
section 2. We emphasize however that the Ndlll lines are 
formed very high in the atmosphere - T5000 ~ 10~^ — 10~^ 
- far above the continuum-forming region tsooo ~ 1 that we 
observe in photometry, hence the two sets of observations 
do not sample the same layers of the star. The Nd ill line 
formation is above the outer turning point of the mode re- 
flection and therefore radial velocity measurements of this 
line do not correspond to the trapped magneto-acoustic 
modes but rather to the running part of mode in the atmo- 
sphere. 

As mentioned in section 2, our formalism is applicable to 
the modes that are weakly distorted by magnetic field such 
that they can be represented by a single value oi £ = 1 . We 
therefore assume that it is the case for HR 3831. Our for- 
malism enables us to determine the mode inclination with 
respect to the rotation axis by comparing with observed 
light curves. One needs an independent determination of 
the magnetic field axis orientation to conclude about the 
alignment of the mode with magnetic axis. 

The magnetic field of HR3831 was first measured by 
Thompson (1981), who found a dipole field and determined 
the star's magnetic ephemcris. These magnetic data and a 
large set of pulsation data were analysed by Kurtz et al. 
(1992), who found that the time of the magnetic minimum 
coincides with the time of the pulsation maximum. More 
precisely, the weakest of the two maxima of the rotational 
envelope of pulsation coincides with the positive extrcmum 
of the magnetic field variations. 

More recent magnetic observations were carried by 
Mathys (1995), Bagnulo et al. (1999) and Kochukhov et 
al. (2004). They found drastically different results concern- 
ing both the polar magnetic strength (Bp) and orientation 
toward the rotation axis (/3). In the first two papers, the val- 
ues of the field strength are considerably stronger {Bp ~ 11 
kG) than that of Kochukhov et al. {Bp ~ 2.5 kG). The 
orientations with respect to the rotation axis are also very 
different with /3 = 7° in Bagnulo et al. (1999) and 87° in 
Kochukhov ct al. (2004). These differences are discussed 
in Kochukhov et al. (2004) who argued that these discrep- 
ancies are probably caused by the neglect of Mathys and 
Bagnulo et al. of the presence of stellar spots in their anal- 
ysis. 

We emphasize at this stage that the magnetic measure- 
ments of Bagnulo et al. (1999) appear incompatible with 
the presence of dipole-like pulsation aligned with the mag- 
netic axis. That would indeed imply an improbably high 
pulsation amplitude that is not believable: The observed 
amplitude of the pulsation would be R{X = 82°) = R{X = 
0°)cos(82). The observed semi- amplitude for IIR3831 is 
about 5 mmag in B; with this geometry, if the mode could 
be seen pole-on (A = 0), then it would have an intrinsic 
semi- amplitude of 36 mmag! This is vastly larger than the 
largest amplitude seen for any of the known roAp star of 
8 mmag, i.e. in the case of HD 60435, and that is for multi- 
ple modes beating with each other. 
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Fig. 4. Applications to four roAp stars: HR3831 (upper left), HD6532 (upper right), aCir (lower left), and HD 99563 
(lower right). The dots represent the observed values and the full lines the theory. Observational data are from Kurtz et 
al. (1997a) for HR3831, Kurtz et al. (1996a) for HD6532, Bruntt et al. (2009) for aCir, and Handler et al. (2006) for 
HD 99563. In our linear treatment of pulsation, the absolute amplitude of the mode is arbitrary. Therefore, the amplitudes 
of the light curves are scaled so that the observed ones equal the theoretical ones at $ = 0. The theoretical pulsation 
phase is with respect to an arbitrary zero-point in time, so is shifted to reproduce the observed pulsation phase at rotation 
phase $ = 0. The best-fits to these curves are {i = 68°, 7 = 91°, V = -7°) for HR3831, {i = 40°, 7 = 70°, V = -15°) 
for HD6532 {i ^ 37°, 7 = 20°, t/- = 0.5°) for aCir, and {i = 44°, 7 = 86°, V = -9°) for HD 99563. The different shapes 
of the amplitudes are caused by different inclinations of the mode axes with respect to the line-of-sight of the observer. 
The four modes are almost linearly polarized. The rotational variabilities of the magnetic fields of HD6532, aCir, and 
HD 99563 are not known; such observations will thus test our model for agreement with the angles we have given. 



In their work. Bigot & Dziembowski (2002) determined 
the magnetic configuration {Bp,P,i) of the star by fitting 
the relative amplitudes of the triplet. With these parame- 
ters, they derived the mode inclination. However, this ap- 
proach leads to two independent constraints, and the au- 
thors had to rely on an independent determination of the 
inclination to the line of sight i. At that time, this value 
was available in Bagnulo et al. (1999), who proposed that 
i = 89°. Using that value of i. Bigot & Dziembowski iden- 
tified a mode that is almost linearly polarized along the 
X-axis (ijj is small) but with an axis of symmetry that is 
steeply inclined to the magnetic or rotation axes, 35° and 
42°, respectively. This result was surprising because it sug- 
gested for the time that modes in roAp may not be aligned 
with the magnetic field axis. This result was questioned 
by Kochukhov (2006), who found a mode nearly that is 
aligned with the magnetic axis from radial velocity mea- 
surements. He argued that the geometrical model of Bigot 
& Dziembowski (2002) is too simple to account for the com- 
plexity of pulsation in these objects. 

We however emphasize that the validity of the formal- 
ism in Bigot & Dziembowski (2002) cannot be questioned 
by this argument, the two studies not being based on the 
same observational constraints. A different value of i leads 



to a completely different value of the mode inclination. 
As a matter of fact, if we assume the value proposed by 
Kochukhov (2006), i.e. i = 68°, the geometrical model fits 
the light curve variations very well for both the amplitudes 
and phases, as seen in Fig.|4l The mode is now inclined 
by 7 = 91° with respect to the rotation axis, i.e. almost 
aligned with the magnetic axis (if we assume the value of 
/3 = 87° given by Kochukhov et al. 2004). This mode is 
still found to be almost linearly polarized, i.e. with a small 
non-axisymmetric components, since ~ 7°. The ratio of 
magnetic to centrifugal shifts of frequencies is /i w —0.8. 
As shown in Fig. 2 of Bigot & Dziembowski (2002), even 
if magnetic and centrifugal shifts are comparable {\fj.\ ~ 1), 
the modes can be almost linearly polarized along the mag- 
netic axis if /3 is large. At the time of the positive mag- 
netic extrema, the angle of the pulsation axis with the 
line of sight is a{0) = 23° and half- rotation period later 
a{TT) = 159°. These results are in good agreement with 
those presented in Kochukhov (2006). 



Since the data of Bagnulo et al. (1999) are unlikely, 
we conclude that the pulsation mode in HR3831 is nearly 
aligned with the magnetic axis. 
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5.2. Application to HD6532 

HD 6532 is one of the few roAp stars to show a polarity- 
reversing dipole- like pulsation mode (Kurtz et al. 1996a); 
the others are HR3831 discussed in the last section, 
HD 99563 (Handler et al. 2006, and below), and HD 80316 
(Kurtz at al. 1997b). The star HD 6532 has a short rotation 
period of only 1.944973 d (Kurtz et al. 1996b), hence a rela- 
tively rapid rotational velocity, similar to that of HR3831. 
The geometry of the star is such that one pulsation pole 
is seen for about two-thirds of the rotation period, and 
the other for one-thirds, with a pulsation phase reversal 
at quadrature. With this geometry, one pole is seen at a 
far more favourable aspect than the other, hence at larger 
amplitude. 

Little is known about the magnetic field of HD6532. 
Mathys & Hubrig (1997) did not detect the longitudinal 
field, but found a significant quadratic field of 22 ± 4G. 
Magnetic measurements are not easy to make for HD 6532, 
given its relative faintness {V = 8.4) and rotationally- 
broadened lines. We are thus unable to compare our pulsa- 
tion solution to the light curves of HD 6532 with indepen- 
dent magnetic data, but we can make predictions about the 
geometry of the star that can be tested by future magnetic 
measurements . 

We modelled the pulsation data for HD 6532 taken from 
Kurtz et al. (1996a). The best fit for we have obtained is 
(i = 40°, 7 = 70°, V' = -15°). The mode axis is closer 
to the line-of-sight than in HR3831, i.e. a{0) = 30° and 
a(7r) = 110°. The difference between the two maxima is 
therefore larger for this star than for HR3831; we found 
that R{0)/R{tt) sa 2.5, as seen in Fig. H With these val- 
ues of a, the observer sees the two hemispheres as the star 
rotates, as in HR3831. The important point here is that 
we can match the rotational amplitude and phase varia- 
tions seen in the light curves of HD 6532 with our model 
and make clear predictions about the geometry that can be 
tested against future magnetic measurements. 

5.3. Application to a Cir 

The star a Cir is the brightest of the roAp stars {V ~ 3.2) 
and has basically a single pulsation mode with a frequency 
of 2.442 mHz {P = 6.8 min) (Kurtz et al. 1994; Bruntt et 
al. 2009). Several other modes of consecutive odd-even de- 
grees are present but with smaller amplitudes (Bruntt et 
al. 2009). These authors derived a larger value of the large 
separation = 60 /iHz than the previous estimate (Kurtz 
et al. 1994). This new value of Aiy, combined with a dras- 
tic downward revision of the effective temperature (7400 K) 
thanks to an interferometric radius determination (Bruntt 
et al. 2008), leads for the first time to consistent astero- 
seismic and Hipparcos parallaxes, in contrast to the result 
found by Matthews et al. (1999). There had been a long 
history of attempts to find the rotation period of a Cir un- 
til it was finally determined by an analysis of the pulsation 
data; the rotation period is 4.4790 ± 0.0001 d (Bruntt et al. 
2009). The rotational sidelobes in the amplitude spectrum 
for this star are small, indicating that there is little modu- 
lation of the amplitude, hence this aspect of the pulsation 
mode does not vary much with rotation. 

That is also true for the magnetic field, for which no 
rotational variations have therefore been detected. Mathys 
& Hubrig (1997) detected no longitudinal field, but found 



a quadratic field strength of 7.5 kC As for HD 6532 above, 
we cannot test our model against independent magnetic 
measurements, but the geometry that we find can be tested 
against future, higher precision magnetic measurements. 

We modelled pulsation data for a Cir taken from Bruntt 
et al. (2009); these are white light data taken by the WIRE 
satellite that have a higher precision than the B data of 
Kurtz et al. (1994), but otherwise have the same rotational 
modulation. We used the value oi i = 37° found by Bruntt 
et al. (2008) by combining an interferometric determination 
of the radius with the determination of v sin i by line spec- 
troscopy. Our best-fit is (7 = 20°, tp = 0.5°). In this case, 
the mode axis is still almost linearly polarized along the 
rotation axis (7 = 20°). The angles made by the pulsation 
axis to the line-of-sight, at the time of the two maxima, are 
close, i.e. a(0) = 17° and a(7r) = 57°. Owing to the small 
tilt of the pulsation axis with respect to the rotation one, 
the observer always sees the same hemisphere and there is 
no phase reversal. The phase is therefore nearly constant. 
The ellipticity ?/; is very small for this star with a value of 
0.5°. This mode is very close to an mu = mode. 

5.4. Application to HD 99563 

The star HD 99563 is a singly-periodic roAp star pulsating 
in a distorted dipole mode that shows both pulsation poles 
over the rotational cycle. Photometrically, it is one of the 
larger amplitude roAp stars with a semi-amplitude as high 
as 6mmag in Johnson B (Handler et al. 2006). Its radial 
velocity variations have the highest amplitudes seen for an 
roAp star, reaching 5kms~^ for some Euii lines, as a con- 
sequence of their strong concentration in small abundance 
spots near the pulsation poles (Freyhammer et al. 2009), 
their stratification high in the atmosphere, and the intrin- 
sically high pulsation amplitude of HD 99563 (for an roAp 
star) . The rotation period for HD 99563 was determined 
by Handler et al. (2006) to be P^ot = 2.91179 ± 0.00007d. 
Some variation in both the amplitude and phase with ro- 
tation, plus the strong £ = 3 components found by Handler 
et al. show this star to have a magnetically distorted dipole 
mode. Freyhammer et al. (2009) showed that some rare 
earth elements are overabundant by more than 10^ solar 
values in spots, and that even H shows a nonuniform sur- 
face distribution as a consequence of the rare-earth element 
patches. The result of this is that the observations - even 
in photometry ~ do not sample the surface uniformly. The 
spatial filter caused by the extreme abundance anomalies is 
probably the source of the apparent distortion of the mode 
determined from photometry. 

We modeled the multisite photometric observations of 
Handler et al. (2006) using the rotational inclination i — 
44° determined by those authors. From this, we found that 
7 = 86° and -0 = —9°, as can be seen in Fig. 4. This is the 
poorest fit of our model to the four stars tested, which we 
believe is a result of the strong distortion of the observa- 
tions by the surface spatial filter. An interesting test of this 
would be to study the pulsation of HD 99563 in Johnson V, 
where there is little rotational variation in the mean light 
(see Fig. 4 of Handler et al. 2006). While the pulsation am- 
plitude will be smaller in Johnson V than B, the star is 
of sufficiently high amplitude to ensure that data of good 
signal- to-noisc ratio can be obtained. We expect that such 
observations will show a weaker distortion. Magnetic mea- 
surements of HD 99563 are consistent with the pulsation 
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geometry (see Freyhammer et al. 2009 for a discussion). 
As for HR3831, the mode axis is nearly aligned with the 
magnetic one (/3 = 88°, see Freyhammer et al. 2009). 

5.5. Comparison of the four roAp stars 

On the basis of this small sample, it appears that there 
is no preferential inclination of the pulsation axis to the 
rotation one since 7 = 91°, 70°, 20°, and 86°. The lack 
of magnetic data prevent any general conclusion regarding 
the alignment of the mode and magnetic axes. However, for 
the two roAp stars (HR3831, HD 99563) for which wc have 
the value of magnetic obliquity we found from the light 
curve analysis that the modes are nearly aligned with the 
magnetic axis (7 « (3). The different behaviour of the light 
curves is explained by the different inclinations of their axis 
of symmetry to the linc-of-sight. We note, however, that 
there is a tendency for the modes to have small ellipticities 
of \ip\ = 7°, 15°, 0.5°, and 9°. A new polarity-reversing roAp 
star found by the Kepler mission, KIC 101095926 shows a 
similar small cUipticity of 1-01 = 4° (Kurtz et al., 2011). 
This means that most of the displacement vectors of these 
modes remain very close to the magnetic plane during the 
pulsation cycle. This may be a common property of roAp 
stars that should be investigated by considering a much 
larger sample than the present one. 

6. Conclusions 

We have presented a formalism to explain the aspect of 
the light curves associated with dipole oscillations in roAp 
stars. In this purely geometrical model, all the complexity 
is reduced to three aspect angles, namely the mode and 
observer inclinations with respect to the rotation axis and 
the polarization angle of the mode. This geometrical model 
can reproduce global mode properties and is particularly 
useful for asteroseismic diagnostics. As we have shown, the 
model works well with the four roAp stars considered in 
this paper. Owing to its simplicity, it has limitations and 
cannot reproduce detailed line profile variation or derived 
accurate field strength. This would need a more realistic 
rotation-magneto- hydrodynamic numerical model. 

We have shown that the light curve variations in roAp 
stars are strongly influenced by the magnetic field strength, 
which is parameterized in this study by a single parameter. 
The reason is that an increase in the field strength pushes 
the mode to be either aligned or perpendicular to the mag- 
netic axis. This consequently changes the aspect of the pro- 
jection of the mode onto the line-of-sight of the observer. 
We have also clearly shown that the coincidence between 
the observed pulsation and magnetic maxima is not evi- 
dence of an alignment of their respective axes as had been 
commonly thought. 

We have applied this formalism to four roAp stars and 
in particular the well-studied roAp star HR3831. For this 
star, we revised the conclusion of Bigot & Dziembowski 
(2002), who derived a highly inclined dipole mode with re- 
spect to the magnetic axis. We do not question their for- 
malism, but rather the magnetic data they used. Indeed, 
in either the light curve or triplet analysis, one needs to 
assume an inclination for the observer. The one used by 
Bigot & Dziembowski (2002) was taken from Bagnulo et 
al. (1999) and led to well inclined mode. Using the more 
recent and the more reliable data of Kochukhov (2006), we 



have found a mode that is nearly aligned with the magnetic 
axis, quite close to the original picture of the oblique pul- 
sator model of Kurtz (1982). This means that the magnetic 
field strength is strong (> 1 kG). In this respect, the value 
of 2.5 kG proposed by Kochukhov et al. (2004) agrees with 
the present conclusions. This strength makes the centrifu- 
gal force inefficient to tilt significantly the mode axis far 
from the magnetic axis. 

The geometrical picture of pulsation presented in Bigot 

6 Dziembowski (2002) and in this paper clearly explains 
the photometric data and leads to mode properties that 
are in good agreement with radial velocity measurements, 
at least in the case of IIR3831. In this model, the inclina- 
tion 7 and polarization ip are two global properties of the 
standing mode and therefore do not vary with depth. This 
is the converse of the result that Kochukhov (2006) found 
for radial velocity measurements of Pr and Nd lines. He 
claimed that the geometrical approach is inadequate since 
it docs not account for this variation with depth. We, how- 
ever, emphasize that the variation in the ratio of the peak 
amplitudes of the triplet found by Kochukhov (equivalently 

7 and ip) indicates the presence of a running wave travelling 
in the atmosphere far above the continuum-forming region 
where the fluctuations of luminosity originate. As shown by 
Saio et al. (2010) with a realistic magncto-hydrodynamic 
pulsation model, the phase of the standing mode is con- 
stant in the deep atmosphere and in the interior, whereas 
the running wave in the atmosphere has a changing phase 
with depth. Our formalism applies to the global standing 
mode whose phase is therefore constant, and the parameters 
7 and "0 can be regarded as global geometrical parameters. 

Acknowledgements. L. B. thanks W. Dziembowski for fruitful discus- 
sions at an early stage of this work and the referee, H. Saio, for his 
constructive suggestions. 
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Appendix A: Special types of light curves 

We apply the general formalism developed in this paper to 
special values of the polarization given hy ip = 0, = tt/2, 
and IV'I = 7r/4. 



A. 0.1. Linear polarization ip = 

We consider the important case of modes linearly polarized 
in the plane (B, ft), characterized by a value of = 0. This 
case, which is even idealistic, is not far from reality since it 
is believed that most of the pulsation modes in roAp stars 
arc almost linearly polarized. 

These dipole modes ('0 = 0) always show extrema in 
phase with the magnetic extrema because their pulsation 
axes lie in the plane (B, fl). Their amplitude is written 



i?(<i>) oc |cosa($) 



(A.l) 



To precisely examine the coincidence between magnetic 
and pulsation maxima, we have to consider four angles 
which are the angles between the mode axis and the line- 
of-sight at the two extrema, a($ = 0) and «($ = tt), and 
the angle between the magnetic axis and the linc-of-sight 
at the same phases, A($ = 0) and A($ = tt). 

The conditions required to have maxima at 3> = [tt] , 
i.e. Ea. [3TH32l can be divided into the following three cases: 

— Case I.: the mode shows two maxima at $ = and 
$ = TT, with phases ^' of opposite signs: 

TT 



Case III.: the mode shows one maximum at $ = tt and 
a minimum at <& = 0, with the same signs for the phases 
* : 



|«(7r)| < |a(0)| < - 



(A.6) 



(A.7) 



Cases II and III correspond to orientations of the system 
where the observer sees only one hemisphere as the star ro- 
tates. Examples of cases I and III are illustrated in Fig. lA.ll 

We also note that for some specific inclinations of either 
the rotation axis or the mode axis, the light curve shows 
some extrema of the same amplitude. In the case of i = 
0[7r/4] and/or 7 = 0[7r/2], we have |q:(0)| — |Q:(7r)|, and 
thereby i?(0) = R{tt). 

The coincidence between the magnetic and pulsation 
maxima can occur in many cases depending on the mode 
and the magnetic field orientation. For example, the condi- 
tion to have pulsation maxima at $ = and at $ = tt and a 
magnetic positive maximum at $ = together with a neg- 
ative magnetic minimum at $ = tt combines the relations 
Ea. lA.2l and Eq.l^H which can be written 

\i-l\<^<\i + l\ and < I < N + /3|.(A.8) 

This relation is not very restrictive and allows a wide range 
of values of (i, 7, 



A. 0.2. Linear polarization — 7r/2 

In the case of linearly polarized modes orthogonal to the 
plane (B,r2), i.e. with ~ tt/2, wc have 



R{^) oc I sin$| 



(A.9) 



No possible coincidence between magnetic and pulsation 
maxima can occur because the mode pulsates always in a 
direction orthogonal to the magnetic axis. The magnetic 
and pulsation maxima are always in quadrature. 

A. 0.3. Circular polarization IV^I = n/A 

In the case of m = ±1 modes in the rotation reference 
system, (ip = ±7r/4, 7 — 7r/2), the light curve simplifies to 
the simple expression 



|a(0)| < 



< 



(A.2) 



i?($) = |sini| ±$. 



(A.IO) 



Hn)\ <-< |a(0)|. 



(A.3) 



In this case, the observer sees the two hemispheres of 
the dipole as the star rotates. The light curves resemble 
those shown in Fig. lA.ll 

Case II.: the mode shows one maximum at 3> = and a 
minimum at $ = tt, with the same signs for phases ^' : 

|a(0)| < \a{n)\ < ^ (A.4) 



or 



^ < |a(0)| < |a(7r)|. 



(A.5) 



The phase of the light curve then takes the characteristic 
shape of straight lines of slope ±1. This situation appears 
in Fig.H 

Appendix B: Nature of extrema at <I> = <I>o 

The condition to ensure that the maxima are located at 
$ = $0 is written 



(t^ - l)cos2$o < 



(B.l) 



The maxima correspond to the instants when the y-axis of 
the ellipse is the closest to the plane (i, J^)- We emphasize 



12 



Bigot and Kurtz: Theoretical light curves of dipole oscillations in roAp stars 




n 2n 371 4ti ,i 2tt 37t 4tt 

♦ rot = 27TQ: (rod) *rot = 27TQ: (rod) 



Fig. A.l. Representation of the light curves for linear polarized (-0 = 0) dipole modes with different inclinations of the 
pulsation axis. This illustrates the cases discussed in section A. 0.1. In both figures, the position of the observer is fixed, 
i = 40°. Left panel: the dipole axis remains close to the line-of-sight during the rotation of the star. The observer only sees 
one hemisphere of the mode, and thereby does not observe a phase jump during the rotation period. Hence, a maximum 
and minimum of the rotational envelope of the pulsation are observed at <I> = and $ = tt. Right panel: the dipole axis 
is now sufficiently inclined with respect to the line-of-sight so that the two hemispheres of the dipole can been seen by 
the observer. A phase jump by tt radians occurs when the node of the oscillations crosses the plane (i, $7). In this case, 
the amplitude of the rotational envelope leads to two maxima of different amplitudes (q;(0) < a^Tr)). 



here that because of the orientation of the system (rotation 
axis and mode), these extrema are not located exactly at 
$ = 7r/2 [tt], as we might think intuitively, but rather are 
phase shifted, except for 7 = 7r/2 or j = 7r/2 for which they 
are located at $ = 7r/2 [tt]. Their amplitudes relative to the 
extrema at $ = are 

which do not vanish as long as the mode is not linearly 
polarized in the plane (B, fl), i.e. r = 0. 



i?(0) 
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